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Abstract 

Mechanical  control  systems  are  a  very  important  class 
of  nonlinear  control  systems.  They  posses  a  rich  math¬ 
ematical  structure  which  can  be  extremely  important 
for  the  solution  of  various  control  problems.  In  this 
paper,  we  expand  the  applicability  of  design  method¬ 
ologies  developed  for  mechanical  control  systems  by  lo¬ 
cally  rendering  nonlinear  control  systems,  mechanical 
by  a  proper  choice  of  feedback.  In  particular,  we  char¬ 
acterize  control  systems  which  can  be  transformed  to 
Hamiltonian  control  systems  by  a  local  feedback  trans¬ 
formation. 


1  Introduction 

Mechanical  control  systems  are  clearly  a  very  large 
class  of  nonlinear  control  systems.  Their  refined  math¬ 
ematical  structure  of  mechanical  control  systems  can 
be  extremely  useful  and  must  be  exploited  in  their  con¬ 
trol  design.  It  is  therefore  very  natural  that  a  wealth  of 
powerful  design  methodologies  have  been  developed  for 
mechanical  systems.  Examples  include  energy  shaping 
methods  [4,  2,  12],  specialized  controllability  notions 
and  tests  [9,  5],  motion  planning  and  generation  [7], 
among  many  others.  See  also  the  monographs  [11,  14] 
for  several  design  techniques  based  on  the  related  no¬ 
tion  of  passivity. 

In  this  paper,  we  broaden  the  applicability  of  design 
tools  for  mechanical  control  systems  to  other  classes  of 
nonlinear  control  systems  by  proper  choice  of  feedback. 
More  specifically,  we  will  solve  the  following  equiva¬ 
lence  problem:  Given  a  control  system ,  determine  if 
it  is  possible  to  transform  it  to  a  Hamiltonian  control 
system  by  a  feedback  transformation.  We  recall  that 
this  has  been  considered  one  of  the  open  problems  in 
the  area  of  mechanical  feedback  control  systems  as  de¬ 
scribed  in  the  following  passage  from  [3]:  “  Find  other 
techniques  which  enable  one  to  me  feedback  control  for 
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mechanical  or,  indeed,  nonmechanical  systems,  which 
leave  or  put  the  system  into  Hamiltonian  or  Lagrangian 
form.  ” 

After  reviewing  some  notions  of  symplectic  geometry, 
in  Section  2,  we  introduce  a  notion  of  Hamiltonian  con¬ 
trol  systems.  In  Section  3  we  provide  a  simple  test  to 
determine  if  a  given  control  system  can  be  rendered 
Hamiltonian  with  respect  to  a  given  Hamiltonian.  In 
Section  4,  we  determine  sufficient  and  necessary  condi¬ 
tions  for  the  existence  of  any  Hamiltonian  and  a  feed¬ 
back  transformation  rendering  a  control  system  Hamil¬ 
tonian.  It  should  be  noted  that  the  feedback  trans¬ 
formations  are  local  in  nature,  and  with  respect  to  the 
canonical  symplectic  structure.  These  geometric  condi¬ 
tions  are  then  illustrated  with  an  example.  In  Section  5 
we  present  some  topics  for  further  research. 


2  Hamiltonian  Control  Systems 

Several  different  models  of  mechanical  control  systems 
abound  in  the  literature  on  control  of  mechanical  con¬ 
trol  systems.  We  will  adopt  what  we  think  to  be  one 
of  the  simplest  such  models:  Hamiltonian  control  sys¬ 
tems.  To  introduce  it,  we  review  some  elementary  no¬ 
tions  of  symplectic  geometry  [6,  10].  A  symplectic  form 
w  on  a  smooth  manifold  M  is  a  two-form  satisfying  the 
following  properties: 

1.  Nondegeneracy:  w(X,Y)  =  0  for  every  X  €  TM 
implies  Y  =  0, 

2.  Closedness:  doj  =  0. 

where  we  have  denoted  the  exterior  derivative  by  d. 
We  also  use  the  notation  ix'w  to  represent  the  contrac¬ 
tion  of  the  two-form  w  with  the  vector  field  X ,  that 
is  (ixw)(y)  —  w(X,Y)  for  any  Y  6  TM,  The  first 
property  of  the  form  w  is  required  to  obtain  an  isomor¬ 
phism  between  TM  and  T*M  from  w.  This  is  achieved 
by  the  correspondence  1  iyw  6  PM  which  allows 
to  associate  a  unique  Hamiltonian  vector  field  Xjj  with 
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any  smooth  map  H  :  M  — ►  R  through  the  equality: 
i  Xhw  =  dH 

Hamiltonian  vector  fields  conserve  energy  (the  Hamil¬ 
tonian  H )  along  their  flows  <pt.  This  fact  can  be  shown 
by  the  following  computation: 

| rtH  =  mx„H) 

=  rt(dH-XH) 

=  rt(ixHu-xH) 

=  4>*t(uj(XH,XH)) 

=  0 

The  nondegeneracy  condition  on  u>  also  implies  that 
the  dimension  on  M  has  to  be  an  even  number,  see  for 
example  [6].  The  closedness  condition  is  required  to 
ensure  that  the  flow  <pt  of  a  Hamiltonian  vector  field 
Xh  respects  the  symplectic  form,  that  is  —  u), 
as  can  be  seen  again  by  simple  computations  where  we 
denote  by  CxHw  the  Lie  derivative  of  w  along  Xh- 

=  dix„w  +  ix„du;) 

=  &*(ddH  +  i*„dtu) 

=  4>t*{  ixHdw) 

We  now  introduce  the  class  of  Hamiltonian  control  sys¬ 
tems  we  will  use  in  this  paper: 

Definition  2.1  Let  M  be  a  smooth  manifold  equipped 
with  symplectic  form  u>  and  let  U  be  the  input  manifold. 
A  control  affine  system  F  :  M  x  U  — +  TM: 

p 

F  =  X  +  J2YiUi 

i—  1 

is  said  to  be  a  Hamiltonian  control  system  with  Hamil¬ 
tonian  H  if  the  vector  field  X  is  Hamiltonian  with 
Hamiltonian  H. 

Within  the  context  of  Hamiltonian  control  systems  one 
could  also  consider  other  models,  for  example,  one 
could  consider  that  the  vector  fields  Yt  axe  also  Hamil¬ 
tonian  [16]  or  even  the  more  general  class  of  port- 
controlled  Hamiltonian  systems  [17].  However,  we  will 
focus  on  this  simple  model  as  the  techniques  developed 
in  this  paper  extend  to  such  cases  in  a  straightforward 
manner. 


3  Achieving  a  given  Hamiltonian 

We  start  by  determining  if  there  exists  a 
feedback  transformation,  that  is,  a  map1 

lIn  the  current  setting  where  control  systems  are  assumed  to 
be  affine  in  the  inputs  it  is  natural  to  restrict  feedback  transfor- 


a(a:)  +  ff(x)v  :  M  x  U  — *  U  with  f)(x)  invertible 

such  that,  the  feedback  transformed  system: 

X(x)  +  YiWfiiWvj  (3.1) 

i  ij 

is  Hamiltonian  with  a  given  Hamiltonian  H.  Prom  ex¬ 
pression  (3.1)  it  is  clear  that  one  only  needs  to  design 
a(x)  to  change  X,  so  we  will  simply  consider  that  f3(x) 
is  the  identity  map  on  U.  This  question  has  the  follow¬ 
ing  simple  answer: 

Proposition  3.1  Let  F  be  an  affine  control  system  on 
a  smooth  manifold  M  equipped  with  symplectic  from  ui 
and  H  :  M  — *  R  a  smooth  map.  There  exists  a  local 
feedback  transformation  rendering  F  Hamiltonian  with 
Hamiltonian  H  iff: 

d H  -  ixw  €  span{iyjW,  iy2w, ....  iypta}  (3.2) 
or  equivalently: 

(d H  -  i^w)  A  iy,w  A  iy2w  A  ...  A  iypw  =  0 

Proof:  Assume  that  a  such  a  feedback  exists,  then 

the  feedback  transformed  system  satisfies: 

‘*+££=1  M*)v,-(*)w  =  dH 
which  we  rewrite  as: 

p 

dH  -  ixw  =  ^2  ai'iYiUJ 

i=l 

clearly  showing  that  (3.2)  is  satisfied.  Conversely,  as¬ 
sume  that  (3.2)  holds,  then  there  are  locally  defined 
smooth  functions  a,-  :  M  — *  K  such  that: 

p 

dH  —  ixcn  =  ^2  Otiy(w 

i=i 

We  now  define  a  by  the  equalities  a*  =  —a,  which 
define  the  desired  local  feedback.  ■ 

While  the  conditions  in  Proposition  3.1  provide  a  quick 
test  to  determine  if  one  can  transform  a  control  system 
to  a  Hamiltonian  one  with  a  specified  H,  they  are  not 
useful  if  one  wants  to  search  for  a  feedback  transfor¬ 
mation  and  also  a  Hamiltonian.  We  devote  the  next 
section  to  this  problem. 


4  Achieving  any  Hamiltonian 

To  provide  a  solution  for  the  general  case  where  no 
Hamiltonian  is  a  priori  specified  we  will  reshape  condi¬ 
tion  (3.2).  We  start  by  making  the  following  additional 
assumption: 

mations  to  affine  transformations  as  this  ensures  that  the  trans¬ 
formed  system  is  still  input  affine. 


1.  The  distribution  spanned  by  the  input  vector 
fields  YuYi,...,Ypt  denoted  by  A,  is  locally  of 
constant  rank. 

4,1  A  Geometric  Solution 

To  develop  a  geometric  solution  we  introduce  the  sym- 
plectic  orthogonal  of  A,  which  we  denote  by  A“  and 
define  by: 

&u  =  {Z£TM  :  w(Z,Y)  =  0  VYeA} 

Note  that  Aw  is  locally  of  constant  rank  in  virtue  of 
nondegeneracy  assumption  and  smoothness  of  u>.  By 
making  use  of  A"  we  can  reformulate  condition  (3.2) 
to  a  more  useful  form.  If  condition  (3.2)  is  satisfied 
then  d H  —  ixw  =  ai'1YlUJ  and  contracting  this 

expression  with  any  vector  field  in  Aw  we  have: 

d H(Z)  =  co(X,  Z)  VZ  £  Aw'  (4.1) 

Conversely,  if  (4.1)  is  satisfied,  then  d H  =  i +  7  for 
some  7  €  spanfi^w,  iy2w, . . . ,  iypw}  which  implies  con¬ 
dition  (3.2)  and  shows  how  (3.2)  can  be  equivalently  ex¬ 
pressed  as  (4.1).  This  expression  can  also  be  regarded 
as  a  partial  differential  equation  (PDE)  whose  solution 
provides  the  desired  Hamiltonian  H.  We  now  interpret 
this  PDE  geometrically  by  defining  the  new  manifold: 

M  ~  M  xR 

with  local  coordinates  (x,  y),  where  x  are  coordinates 
for  M  and  y  coordinates  for  R.  We  also  define  the  map 
H  :  M  ->  R  by: 

H(x,y)  =  H(x)  -  y 
as  well  as  the  vector  fields: 

z,zl  +  .ix,z)± 

defining  the  distribution  A“  on  M.  These  new  objects 
allow  to  rewrite  (4.1)  as: 

CzH  =  0  VZ  €  A^ 

which  we  interpret  as  the  requirement  that  vector  fields 
Z  are  tangent  to  the  submanifold  (H)~1( 0)  of  M. 
We  thus  see  that,  in  this  geometric  interpretation  of 
PDE  (4.1),  finding  a  Hamiltonian  H  is  equivalent  to 
finding  a  submanifold  of  M,  implicitly  defined  by  a 
map  H  :  M  — »  R,  such  that  A"  is  contained  in  the 
tangent  space  of  (,H’)~1(0).  Furthermore  H  must  also 
satisfy  ^  ff  0  which  ensures,  via  the  implicit  function 
theorem,  that  H  defines  a  function  on  M,  the  desired 
Hamiltonian  H.  Necessary  and  sufficient  conditions  for 
the  local  existence  of  such  a  map  H  are  given  in  the 
next  result  which  provides  geometric  conditions  for  the 
integrability  of  PDE  (4.1): 


Theorem  4.1  Let  F  be  an  affine  control  system  on  a 
smooth  manifold  M  with  symplectic  form  oj  and  denote 
by  C  and  C  the  involutive  closures  of  A"  and  A“  which 
we  assume  to  be  regular.  There  exists  a  locally  defined 
map  ff  :  M  -*  1  and  a  local  feedback  transformation 
rendering  F  Hamiltonian  with  Hamiltonian  H  iff: 

dim(C)  =  dim(C)  (4.2) 

Proof:  In  view  of  the  discussion  preceding  Theo¬ 

rem  4.1  it  suffices  to  show  that  (4.2)  is  sufficient  and 
necessary  for  the  existence  of  the  map  H  :  M  — *  R 
satisfying  ^  ff  0  and  Aw  C  T((Pf)‘“1( 0)).  To  show 
necessity  assume  the  existence  of  H  =  H  —  y  satisfy¬ 
ing  Aw  C  T((H)~1( 0)).  Let  i  :  M  — >  M  be  the  map 
i(x)  =  (x,H(x))  and  note  that: 

TJ-Z(x)  =  Z{x)^+TXH -Z{x)~ 

=  Z(x,H(x)) 

=■  Z  o  j(x) 

where  the  second  equality  follows  from  (4.1).  This 
shows  that  the  vector  fields  Z  are  i-related  to  the  vec¬ 
tor  fields  Since  if  Zj  is  f-related  to  Z,  and  Z;-  is  i- 
related  to  Zj,  the  bracket  [Zj,  Zj]  is  i-related  to  [Zj,  Zj] 

( see  [1]),  it  follows  by  induction  that  dim(C)  =  dim(C). 

Sufficiency  is  proved  by  applying  Frobenius  theorem  to 
C  (which  is  regular  by  assumption)  to  conclude  the  exis¬ 
tence  of  a  submanifold  of  M  to  which  the  vector  fields 
in  Aw  are  tangent.  Furthermore,  Frobenius  theorem 
also  ensures  that  this  submanifold  is  locally  described 
by  the  zero  level  of  a  smooth  map  H  :  M  —>  R  [1].  It 
remains  to  show  that  ^  ff  0.  We  proceed  by  contra¬ 
diction  assuming  that  dim(C)  =  dim(C)  and  ^  =  0. 
Then,  the  vector  field  Z  =  Ojy  +  1  €  ker(di?)  = 

kcr(^dx  +  ^d y).  This  shows  that  dim(C)  is  at  least 
greater  them  dim(C)  by  one,  a  contradiction.  ■ 

Theorem  4.1  gives  necessary  and  sufficient  conditions 
for  the  existence  of  a  solution  to  PDE  (4.1).  The  so¬ 
lution  of  partial  differential  equations  naturally  appear 
in  similar  local  feedback  equivalence  problems,  such  as 
feedback  linearization  [8],  and  several  control  design 
problems  for  mechanical  control  systems  [2,  12). 

As  an  immediate  consequence  of  Theorem  4.1  we  see 
that  in  the  case  dim(A")  =  1,  that  is,  control  system  F 
has  2to  —  1  inputs  and  dim(M)  =  2 to,  condition  (4.2) 
is  automatically  satisfied. 

Corollary  4,2  Let  F  be  an  affine  control  system  on  a 
smooth  manifold  M  of  dimension  2m  with  symplectic 


form  to  and  2m  —  1  inputs.  Then,  there  exists  a  locally 
defined  map  H  :  M  — *  R  and  a  local  feedback  trans¬ 
formation  rendering  F  Hamiltonian  with  Hamiltonian 
H. 

The  discussion  so  far  has  assumed  that  the  symplectic 
form  a)  has  been  a  priori  specified.  However,  this  is 
not  necessary  in  virtue  of  Darboux  theorem  [6]  which 
asserts  that,  locally,  every  symplectic  manifold  of  di¬ 
mension  2m  is  sympleetomorphic  (diffeomorphic  by  a 
diffeomorphism  that  preserves  the  symplectic  forms)  to 
M2m  with  symplectic  form: 

m 

to  =  ^  daq  A  dj h  (4.3) 

<= l 

expressed  in  coordinates  (xi,x2,  • .  ■,xm,V’i,y2,  ■  ■  ■  ,Vm) 
for  R2m,  We  thus  see  that  if  a  control  system  can  be  lo¬ 
cally  rendered  Hamiltonian  with  respect  to  the  form  w, 
then  by  a  change  of  coordinates,  it  is  also  Hamiltonian 
with  respect  to  any  other  symplectic  form. 

4.2  An  alternative  Characterization 
The  conditions  for  the  existence  of  a  Hamiltonian  and  a 
feedback  transformation  given  in  Theorem  4.1  require 
the  computation  of  several  objects  such  as  Au,  A “,  C, 
C,  etc.  However,  some  of  these  objects  contain  some 
degree  of  redundancy  and  we  will  now  see  how  one  can 
verify  the  conditions  of  Theorem  4.1  in  a  more  direct 
way.  In  particular  we  shall  take  advantage  of  the  special 
form  of  the  vector  fields  in  A“.  From  the  expression  of 
[Zi,Zj\  in  local  coordinates: 

=  [d u(X,Zj)  Oj  Lo{X,Zi) 

_  '  1 1  o  f 

d to(X,Zi)  oj  to(X,Zj) 
we  see  that  [Zi,Zj]  is  given  by: 

[Zi,  Zj}—  +  (. CzMX ,  Zj)  -  CZjoj{X,  ZCi)  ~  (4-4) 

We  now  rewrite  £ztco(X,  Zj)  —  £Zjlo(X,  Zfi)  as: 

-d{izfixaj)(Zi)  -  CZju{X,  Zf) 
which  by  the  Cartan  magic  formula  [1]  becomes: 

-(£z3- ixto  -  izj<iixto)(Zi)  —  £Zjio(X,  Zf) 

—  -izffizffxto  +  iziizjdiyw  -  LZjto{X,  Zfi) 

—  Z-zfizfixw  -  iziC-zfixa)  +  izfizjdixu 
=  w(X,  {ZuZjV+izfiz.dixw  (4.5) 

and  where  the  last  equality  is  a  consequence  of  the  fact 
that  CzfiZil  -  i =  7 ([Zi,  Zj})  for  any  one-form 
7- 

Expression  (4.5)  allows  to  formulate  the  following  al¬ 
ternative  version  of  Theorem  4.1: 


Theorem  4.3  (Alternative  Characterization) 

Let  F  be  an  affine  control  system  on  a  smooth  mani¬ 
fold  M  with  symplectic  form  to.  There  exists  a  locally 
defined  map  H  :  M  — *  R  and  a  local  feedback  transfor¬ 
mation  rendering  F  Hamiltonian  with  Hamiltonian  H 

iff: 

izj izj dixto  =  0  VZ{,  Zj  €  A"  (4.6) 

Proof:  As  we  have  seen  in  the  proof  of  Theorem  4.1, 
existence  of  H  and  the  feedback  transformation  implies 
that  every  vector  field  Z  6  A"  is  i-related  to  the  vec¬ 
tor  field  Z  €  Aw  for  i(x)  =  ( x ,  H{x)).  This,  in  turn, 
implies  that  [Zi,  Zj]  is  also  i-related  to  [Z, ,  Zj]  leading 
to: 

=  [Zi.Zj]  ^  +»(X,  [ Zi,Zj ])A 

Comparing  this  expression  with  (4.4)  we  see  that  we 
necessarily  have  iZiiZjdixco  =  0. 

To  show  sufficiency  we  note  that  if  [Zi,  Zj]  = 

for  Zi,  Zj,  Zi  €  C  and  smooth  real  valued  functions  o;, 

it  follows  by  (4.4)  that  [Zi,  Zfi]  equals: 

[Zi,Zj}~  +  {uiiXfiZi'Zj})  +iZiiZjdixw)^ 

=  YYa^ffx  +  HX^Z)  +  izfi^dixto)^ 

Q  Q 

=  YL aiZt  &  +  ( Y1  aiuj(x’ Zi ) + iz<  di*  w)  gf; 

=  Yj  a‘  Zl  +  ‘Zi  lzi  di  Jf 1 w 

The  assumption  izfiZjdixto  —  0  now  allows  to  con¬ 
clude  that  \Zi,Zi]  =  atZi  =>  [Zi,Z[]  =  'ffio-iZi 

and  an  induction  argument  shows  that  dim(C)  = 
dim(C)  which  by  Theorem  4.1  implies  the  existence  of 
H  and  the  desired  feedback  transformation.  ■ 

This  formulation  also  allows  to  see  that  Corollary  4.2 
is  a  simple  consequence  of  the  skew-symmetry  of  w. 
If  dim(A“)  =  1  we  have  that  for  any  Zi,  Zj  €  A“, 
Zi  =  A  Zj  for  a  smooth  real  valued  map  A  and 
(dixuj)(Zi,  XZj)  =  X(dixto)(Zi,  Zfi)  =  0. 

The  necessary  steps  to  determine  the  existence  of  a 
solution  to  PDE  (4.1)  can  now  be  resumed  to  the  fol¬ 
lowing: 

1.  Compute  a  basis  for  the  sym¬ 
plectic  orthogonal  Aw  of  A, 
A “  =  {Z£TM  ;  u[Z,Y)  =  0  VY  €  A}, 

2.  Check  if  iZiiZjdixto  —  0  for  every  Zi,  Zj  in  the 
basis  of  A"  computed  in  step  1. 


which  by  differentiation  gives: 


If  Theorem  4,1  is  satisfied,  then  a  solution  to  (4.1) 
must  be  obtained  in  order  to  determine  the  feedback 
transformation.  If  H  is  such  a  solution,  we  determine 
the  feedback  transformation  by  computing  the  smooth 
functions  ai  satisfying: 

p 

d  H  —  ijfta  — 

i=  i 

These  functions  allow  to  determine  the  term  a(x)  of  the 
feedback  transformation  a(x)  +  / 3(x)v  by  the  equalities 
a,  =  —a,.  The  term  Pix)  can  be  taken  as  the  identity 
on  U  or  any  other  invertible  (pointwise)  linear  map 
from  U  to  U. 

4.3  Example 

We  now  provide  an  example  of  the  previously  intro¬ 
duced  methodology.  Consider  the  following  control  sys¬ 
tem: 

±i  =  X2X3+X2U! 

x2  =  x\x\  +  -x\ 

x3  =  X1X2X4  +  X2X3U2 
x 4  =  xix3  (4.7) 


di x<jJ  =  —  2xix2dx2  A  dxi  —  x3dx3  A  dxi 

+x2dx3  A  dx2  —  x3dxi  A  dx3 
+x2x4dx!  A  dx4  +  X!X4dx2  A  dx4 
=  —  2xix2dx2  A  dxi  +  x2dx3  A  dx2 

+x2x4dx!  A  dx4  +  X1X4 dx2  A  dx4 

and  evaluating  dl xu  on  Z\  and  Z2: 

iz2dixoJ  =  x2dx2 
iziiz2dixw  =  0 

shows,  via  Theorem  4.3,  that  a  Hamiltonian  and 
a  feedback  transformation  exist. 


To  obtain  H  one  has  to  solve  (4.1),  which  in  this  case 
results  in: 


dH 

dxi 

dH 

dx3 


1  2  2 

~2X3  -  XlX2 

-X1X3 


It  suffices  to  solve  the  first  equation  to  obtain: 


on  K4  with  symplectic  form  w  =  dx\  Adx2+dx3  Adx4. 
In  this  case  we  have: 

X2X3  x2  0 

X  _  XXX\  +  \x\  y  _  0  y  _  0 

xjx2x4  1  0  2  x2x3 

X1X3  0  0 

(4.8) 


H  =  ~^(xjx%  +2:1X3) 

One  now  computes  dH  —  i xw  =  (xfx2  -  x2X3)dx2  — 
xjx2x4dx4  which  can  be  written  as  dH  —  ixoj  = 
aiiyjtu  +  a2iy2w  for: 

ai  =  (x3-x?)  a2  =  —  — X1X4  (4-11) 

2:3 


equation,  and  provided  sufficient  and  necessary  condi¬ 
tions  for  the  local  existence  of  solutions.  These  results 
enlarge  the  class  of  systems  to  which  powerful  control 
design  methods  developed  for  mechanical  systems  are 
applicable. 

Many  related  problems  remain  open.  When  we  can¬ 
not  perform  such  a  feedback  transformation  it  may  still 
be  possible  to  extract  a  quotient  (an  abstraction,  see 
for  example  [13,  15])  or  a  subsystem  that  is  mechan¬ 
ical,  or  that  can  be  rendered  mechanical  by  feedback. 
This  would  allow  to  synthesize  controllers  for  part  of 
the  variables  by  making  use  of  techniques  developed 
for  mechanical  control  systems.  Another  alternative  is 
to  perform  a  change  of  coordinates  on  the  state  space 
that  would  render  the  system  Hamiltonian,  possibly 
complemented  by  a  feedback  transformation. 
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